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Tendex and vortex fields, defined by the eigenvectors and eigenvalues of the electric and magnetic 
parts of the Weyl curvature tensor, form the basis of a recently developed approach to visualiz- 
ing spacetime curvature. In analogy to electric and magnetic fields, these fields are coordinate- 
dependent. However, in a further analogy, we can form invariants from the tendex and vortex fields 
that are invariant under coordinate transformations, just as certain combinations of the electric and 
magnetic fields are invariant under coordinate transformations. We derive these invariants, and 
provide a simple, analytical demonstration for non-spherically symmetric slices of a Schwarzschild 
spacetime. 
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Tendex and vortex fields have recently been introduced 
as tools for the visualization and interpretation of gravi- 
tational fields and spacetime dynamics (pQ, see also [2]). 
In particular, it has been suggested that they can pro- 
vide deeper insight into phenomena observed in numeri- 
cal simulations of binary black hole mergers, including or- 
bital hang- up (see, e.g., [3]), spin-flip (e.g. [4]) and black 
hole recoil (e.g. [SHS]). Tendex fields are defined in terms 
of the eigenvalues and eigenvectors of the electric part of 
the Weyl curvature tensor and describe tidal stretching 
or compression, while vortex fields are defined in terms 
of the eigenvalues and eigenvectors of the magnetic part 
of the Weyl tensor and describe precession. Examples of 
tendex and vortex fields for a number of different space- 
times can be found in references [TJ H0Hl2] . 

In analogy to their electromagnetic counterparts, ten- 
dex and vortex fields are not invariant under general co- 
ordinate transformations. However, just as for electro- 
magnetic fields, it is possible to find combinations of the 
tendex and vortex fields that are invariant. In this short 
paper we derive these invariants and also provide a simple 
example that illustrates both the coordinate dependence 
of the tendex and vortex fields themselves, as well as the 
coordinate-independence of these invariants. 

Before considering tendex and vortex fields it is useful 
to review the properties of electromagnetic fields under 
coordinate transformations. Recall that the electric field 
E l and the magnetic field B l are given by components of 
the four-dimensional, antisymmetric Faraday tensor F ab 
(here and in the following, indices i, j, etc. denote spatial 
indices, while indices a, 6, etc. denote spacetime indices). 
Under a general coordinate transformation, the Faraday 
tensor transforms as a rank-2 tensor, and the new electric 
and magnetic fields can be identified from the Faraday 
tensor in the new coordinate system. In particular, a 
coordinate transformation may therefore mix the electric 
and magnetic fields. 

More specifically, we may write the Faraday tensor as 

pab = n a E b _ n b E a + ^dabc^ ^ 

Here we have assumed a foliation of the spacetime M into 
a family of spatial slices X, each of which corresponds, 



at least locally, to a level surface of a coordinate time t 
(see, e.g., [13]). The normal vector n a is defined as the 
future-pointing normal on the spatial slices, and e a bcd is 
the Levi-Civita tensor. A normal observer then identifies 

E a = F ab n b , B a = ~e abcd F cd n b . (2) 

as the electric and magnetic fields. By construction both 
vectors are spatial, i.e. E a n a — B a n a = 0. 

Under a general coordinate transformation, level sur- 
faces of the coordinate time t may change, meaning that 
the new spatial slices may be tilted with respect to the 
old spatial slices (as is familiar from boosts in special rela- 
tivity) . This means that the normal vector n a in the new 
coordinate system does not represent the same object as 
that in the old coordinate system - it points in a different 
direction. From the relations ([2| it is therefore clear that 
E a and B a do not transform as four-dimensional rank-1 
tensors under general coordinate transformations (even 
though they do transform as three-dimensional vectors 
under purely spatial coordinate transformations, which 
leave n a unaffected). 

However, the electric and magnetic fields inherit in- 
variants from the Faraday tensor. In particular, from 
the invariance of the scalar F ab F ab it follows that 

I™ = E 2 - B 2 (3) 

must be invariant. Here E and B denote the magnitudes 
of the electric and magnetic fields. A second invariant 

lf M = E a B a = E,B l = EB cos 6 EB (4) 

can be derived from the invariance of the determinant of 
the Faraday tensor (under transformations between or- 
thonormal frames). Here 9eb denotes the angle between 
E 1 and B l . It can also be shown that any other invariants 
for the electromagnetic fields are combinations of these 
two (see, e.g., Problem 4.2 in [14]). 

In general relativity, we may similarly write the Weyl 
curvature tensor C a bcd in terms of its electric part £ ab 
and its magnetic part B ab , 

Cabcd — 4 (ri[ a n[ c + 7[ [c) £&]d] + 2e abe n[ c B d ] e 

+ 2n [a B b]e e cd e , (5) 
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where e a bc = n d edabc is the spatial Levi-Civita tensor, and 
where jab is the spatial metric induced on the spatial slice 
(see [H]). A normal observer then identifies the electric 
and magnetic parts as 



£ a b — C ac bdn c n d , B a b — — *C, 



acbdn n 



(6) 



where *C a bcd = \tab e ^ C e f ca < denotes the dual of the Weyl 
tensor (compare PS]). By construction, both £ a b and B a b 
are spatial, symmetric and traceless. By the same argu- 
ments as above we also note that £ a b and B a b do not 
transform as four-dimensional rank-2 tensors. Their be- 
havior under general coordinate transformations instead 
follows from that of the rank-4 Weyl tensor C a bcd- 

In analogy to their electromagnetic counterparts, £ a b 
and B a b inherit invariants from the Weyl tensor. In vac- 
uum, the Weyl tensor has four independent algebraic real 
invariants, which may be written as the real and imagi- 
nary parts of the scalars 



Jf, [CabcdP 



cdab 



;/-~r * s~icdab\ 
l^abcd l> J 



(7) 



and 



Jtri flcdef si ab • St /-icdef '* /~i ab\ /o\ 
= QQ\ abcd e f ~ lL/ abcd<^ O e / ), (O) 

( [H] 1 see also [HI US] ) ■ These well-known scalars play a 
central role in the Petrov classification of spacetimes (see, 
e.g., [Il])- They also form the basis of some spacetime 
diagnostics, including the specialty index suggested by 
|2"0] , and a radiation scalar suggested by [TB] ■ In [T^] , they 
are used to identify a geometrically motivated coordinate 
system for numerical-relativity simulations. In terms of 
Sij and Bij , the invariants I and J can be written as 



- (£ £ lJ 
2 \°v c ' 



(9) 



and 



J 



- l -£ l P k £\ + ^£ i jB j k B k i 
i (^B i j B j k B k i - ^ s ^ k e\ 



(10) 



(see EI], US]). 

We now express these invariants in terms of tendex and 
vortex fields. Tendex fields are defined by the eigenvalues 
and eigenvectors of the electric part of the Weyl tensor, 
£ij, while vortex fields are defined by those of the mag- 
netic part Bij . The eigenvalues A and eigenvectors v % are 
most easily found in an orthonormal basis (denoted with 
hats) 



(11) 



and similar for Bij. A negative (positive) tendex eigen- 
value indicates that observers will tend to be tidally 
stretched (compressed) in a direction aligned with the 



corresponding tendex eigenvector, while a positive (neg- 
ative) vortex eigenvalue indicates that a gyroscope will 
exhibit (counter) clockwise differential precession relative 
to a nearby gyroscope in the direction of the correspond- 
ing vortex eigenvector [TJ [TU] . 

Since both £ a b and B a b are spatial and symmetric, they 
each have a set of three orthonormal eigenvectors. We 
may therefore write the two tensors in terms of their 
eigenvalues A and orthonormal eigenvectors v l as 



£ l ° = \ E x{v Ex y{v Ei y + XE2{VE2T{VE1) 3 



^ XEk{vEk) l {v E k) J 



k=l 



and similarly 



b^ = Y, x Bk(vBkT(v Bk y. 



(12) 



(13) 



k=l 



We use explicit summations for sums over eigenvalues 
and eigenvectors, but will continue to use the implied 
Einstein summation rule for sums over vector indices. 

We can now derive invariants for the tendex and vortex 
fields by inserting these expressions into the invariants I 
(eq. (pj)) and J (eq. ( 10 ) ) . From the real part of 2" we 
find the first invariant 



1 3 

^1 = o /, ( X Ei - X Bi) ' 



(14) 



in close analogy to the electromagnetic invariant pj. 
Likewise, the second invariant follows from the imaginary 
part of I, 



i=l j=l 



3 3 



/ t / , ^Ei^Bj COS Q E iBj j 
i=l j = l 



(15) 



and closely resembles the electromagnetic invariant Q. 
The last two invariants, 



Ji 



1 3 

i=l 



^33 

+ X Ei^Bj COS 2 e BiB j (16) 



and 



1 3 

j% = g zl x b* 

i—l 



i=l j = l 



-^3 3 

n /] X Bi^Ej COS 2 9 E iBj , 
i=l j=l 



(17) 



can be derived from the real and imaginary parts of J . 

We note in passing that only three of the nine angles 
0EiBj between (vEi) k and [y B j) k are independent. It is 
possible to write these nine angles in terms of three Eulcr 
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angles, but that does not appear particularly useful in 
this context. 

As a simple analytical demonstration we now con- 
sider a Schwarzschild spacetime. On a slice of constant 
Schwarzschild time t, the spatial metric, expressed in 
Schwarzschild coordinates, is 



jij = diag ((1 - 2M/R) 1 ,R 2 ,R 2 sin 2 



(18) 



while the extrinsic curvature Ay vanishes. In a spheri- 
cal polar orthonormal basis the electric part of the Weyl 
tensor is 



% = diag (-2M/R 3 , M/R 3 ,M/R 3 



(19) 



(see, e.g., [13 [22]), so that we can immediately identify 
the eigenvalues 



A 



El 



2M M 
and can choose the corresponding eigenvectors to be 



(20) 



( VE1 y = {e k )\ (v E2 y = ( e§ )\ ( VE3 y = ( e<j> y. (21) 

The magnetic part of the Weyl tensor By and its associ- 
ated vortex fields, meanwhile, vanish identically on slices 
of constant Schwarzschild time. 

Now consider a different slicing of the same spacetime. 
The new slices are level surfaces of a new time coordinate 
t, which we constructed from the old time coordinate t 
with the help of a "height function" h, 



t = t + h. 



(22) 



We will introduce a coordinate system that breaks spher- 
ical symmetry by choosing ft, to be a function of 8, 
h = h(9). Transforming to the new coordinate system, 
we can now identify the spatial metric and extrinsic cur- 
vature on a slice of constant time t (see [12] for a peda- 
gogical example). The spatial metric is given by equation 
(18), but with jgg = R 2 replaced with 



lee 



R\ 



where we have defined 



77 = 1- 



{dehf 
R 2 



1 - 



2M\ 
~R~) 



(23) 



(24) 



When h is constant, we have r\ = 1 and recover results 
for slices of constant time t. The extrinsic curvature no 
longer vanishes; its nonzero components are now 



A" 



ne 



= K eR = 



Kee = 



,1/2 



1 (d e h) 3M-R 
V 1/2 R 2 ^1-2M/R' 

(d 2 h) y/l - 2M/R, 



and 



Ka,a, = 



,1/2 



(d e h) y/l - 2M/R 



sin tf cos ( 



(25) 
(26) 

(27) 



In a spherical polar orthonormal basis the electric part 



of the Weyl tensor is given by equation (19), but with 
S~~ = —2M/R 3 replaced with 



RR ~ R 3 { r, 
and with Eii = M/R 3 replaced with 



£ xx 



M 
R 3 



V 



(28) 



(29) 



Since £y is still diagonal, the eigenvectors are still given 
by equations (21 ), and the corresponding eigenvalues are 



Abi - £ RR , 



Afi2 - £§§ - -^3 , A_B 3 - S^, (30) 



with £ RR and £^ given by equations ( 28 ) and ( 29 1 , re- 
spectively. 

For a slice of constant time t, the only non-vanishing 
components of the magnetic part of the Weyl tensor are 



a R<t> °4>R £>3 



The eigenvectors of By are 



and 

(v B3 y 

with eigenvalues 

Asi = — 

and 



1 



V 



{e R f + (ei)- 



(W = (e § y 



VB3 



3M 



A 



B2 



0. 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



This example highlights the fact that tendex and vortex 
fields are slicing-dependent concepts. 



We are now ready to calculate the invariants ( 14 )-( 17 > 



in both coordinate systems. It is easiest to start with a 
slice of constant Schwarzschild time t. In that case, By 
vanishes, so its eigenvalues and eigenvectors do as well. 
Substituting equations (20) into equations (14)-(17), we 
readily find 



3M 2 



Jl 



M 3 



J 2 = o. 



(37) 



On the slice of constant coordinate time t we have 
A732 = 0, so that the invariant I 2 reduces to 

3 

1-2 = ^2 XEi ( Asi cos2 9eiB1 + Xb3 cos2 °EiBa) ■ (38) 
i=l 
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But we have cos 2 6 E iB3 — cos 2 0EiBi for each eigenvector 
(vEi) k , as well as Xb3 = — Xbi, so that I2 again vanishes. 
The argument for J2 is very similar. 

To evaluate the invariant X\ on a i = const slice we 
first observe that 



2 _ 6M 2 

i—1 



ism 2 1 

~R^ 



V 



(39) 



Adding the contributions of the eigenvalues \ 2 Bi from 
equations ( 35 ) and ( 36 1 we see that the second term in the 
above expression cancels exactly, leaving X\ = 3M 2 /R 6 
as before. 



To calculate J7i, we note that the first term of (16) 
evaluates to 



1^,3 _M 3 

A Ei — 



6^ 

■I— 1 



R 9 



9(1-7?) 
2rf 



(40) 



Carrying out the sum over j in the second term of (16 1, 
and using Xb2 — 0, we find 



After considering the angles between the eigenvectors and 
using the fact that the electric part of the Weyl tensor is 
traceless, we see that (41) can be written as 



iff) A 2 , M3 9(1-7?) 

2 2^ 2^ Xsi X Bj co& OziBj - - -jp 2 2 • 

i=l 3=1 ' 



(42) 



Adding this to equation ([40]) shows that the invariant 
Ji remains J\ — M 3 /R 9 , which is not surprising but 
reassuring. 

In conclusion, we have derived slicing-invariant quan- 
tities for tendex and vortex fields from well-known scalar 
curvature invariants. These invariants play the same role 
as the invariants for electric and magnetic fields in elec- 
tromagnetism. We also demonstrate both the coordinate- 
dependent nature of the tendex and vortex fields, as well 
as the coordinate-independence of our new invariants, for 
non-spherical slices of a Schwarzschild spacetime. 



^ 3 3 

2 ^ ^ XEiX%j COS 2 0EiBj 
i=l j=l 



9 M 2 1 — 7? 
r? 2 



x x Ei (cos 2 6 Eim + cos 2 Oems) ■ (41) 
i=l 
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